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We calculate the dilepton emission rate from a fireball created in an ultra-relativistic 
■ heavy-ion collision. For the partonic phase, we complement the perturbative results by a 
phenomenological approach based on lattice QCD results, whereas in the hadronic phase 
^} \ temperature and finite baryon density effects on the spectral function are considered. The 
CN ' resulting rates are compared to data from CERES/NA45. 
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q ! 1. Introduction 

From lattice simulations it is known that the theory of strong interactions, QCD, un- 
dergoes a phase transition at a temperature Tq ~ 160 MeV ffl] from a confined hadronic 
phase to a phase where quarks and gluons constitute the relevant degrees of freedom, the 
J^l ■ quark-gluon plasma (QGP). It is hoped that it is possible to create this QGP in ultrarel- 
ativistic heavy-ion collisions at CERN and RHIC. Dileptons (e + e~ and pairs) are 

interesting probes in this context since they do not interact strongly but rather escape 
unthermalized from the hot and dense region formed in such collisions, the fireball. In 
the early stage, a thermalized QGP phase is expected to exist. There, dileptons originate 
mainly from qq annihilation processes. As the fireball expands, it cools off and hadroniza- 
tion sets in at To- In the hadronic phase, the main dilepton sources are pion and kaon 
annihilation processes which are enhanced due to the formation of the light vector mesons 
p, uj and 0. 

The differential dilepton emission rate from a hot domain in thermal equilibrium is given 
by 

dN Q! 2 1 - . . 

-Imr%,T) (1) 
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where II denotes the spin-averaged time-like self energy of a photon in the heat bath, 
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q is the photon four- momentum, a = f-, and the lepton masses have been neglected. 
This rate has to be integrated over the time evolution of the fireball volume V(t) and 
temperature T(t). 
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2. The QGP phase 

Thermal perturbation theory is presumably not applicable to evaluate Imfi in the QGP. 
The coupling constant g s ~ 2 is large, and the maximal temperatures reached in ultra- 
relativistic heavy ion collisions are still 0(Tc), indicating that non-perturbative physics 
is important. Consequently, we try to use a more phenomenological approach to thermal 
QCD. It has been shown that it is possible to describe the equation of state (EoS) of QCD 
as obtained from lattice calculations by the EoS of a free gas of quasiparticles with ther- 
mally generated masses m(T) ||. The model employs 'effective', non-interacting massive 
quarks and gluons as the relevant degrees of freedom, all dynamics is incorporated in a 
T-dependent effective coupling constant g s (T). This coupling is based on the expression 
for the running coupling at T = 0. For the quark masses m q {T) the two parameters T s 
and A are fitted to the lattice data: 

m 2 AT) = -g 2 s (T)T 2 with g 2 s {T) = — (2) 

qK ' r sV ' sK ' (11N C -2N F ) In [{T + T S )/{T C /X)} 2 

Hence, the only QGP contribution to Imfi in eq.([I|) is the gg-loop, evaluated with the 
thermal quark masses (|2]). 

3. The hadronic phase 

Below Tc, confinement sets in and the effective degrees of freedom change. The photon 
couples now to hadronic J p = 1~ states: the p, u and <fi mesons and multi-pion states 
carrying the same quantum numbers. The electromagnetic current-current correlation 
function can be connected to the currents generated by these mesons using an effective 
Lagrangian which approximates the SU(3) flavour sector of QCD at low energies. The 
appropriate model for our purposes is the improved Vector Meson Dominance model 
combined with chiral dynamics of pions and kaons as described in || . Within this model, 
the following relation between the imaginary part of the irreducible photon self-energy 
Imfi and the vector meson self-energies Hv{q) i n vacuum is derived: 

^ lmn^(g) V~%) q2 ~ m2y m 

Imr%) = 2^ 2 F v{q) , F v{q) = ~y 2 , - T rr r \ > 3 

V 9v q 2 - my + tlmU v (q) 

where my are the (renormalized) vector meson masses, g v is the 'yV coupling and g the 
$y coupling, where $ stands for one of the pseudoscalar Goldstone bosons 7i ± ,ii and 
K ± ,K°. Eq.(|3|) is valid for a virtual photon with vanishing three-momentum. Finite 
temperature modifications of the vector meson self-energies are calculated using thermal 
Feynman rules, the explicit calculations for the p-, to- and 0-meson can be found in ||]. At 
SPS energies, there is still considerable stopping of the interpenetrating nuclei, resulting 
in a net baryon density in the central rapidity region. For the evaluations of these finite 
baryon density effects, we use the results discussed in 0. There it was shown that in 
the linear density approximation Ily is related to the vector meson - nucleon scattering 
amplitude T VN 

U v (q°,q = 0;p) = IP^ - P T VN (q), T VN (q) = —J d A x e^{N\Tj,{x)f{Q)\N), (4) 
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with \N) being the nucleoli states with vanishing momentum. In the following, we assume 
that the temperature- and density dependences of Uy factorize, i.e. we replace in eq.(|]) 
Ily ac by the temperature-dependent ITy(T) and leave T VN unaffected. This amounts to 
neglecting contributions from matrix elements like (Nir\Tj fl (x)j ll (0)\TTN) (meson-nucleon- 
pion scatterings where the pion comes from the heat bath). Furthermore, this approx- 
imation does not take into consideration a possible T-dependent pion or nucleon mass. 
The corresponding spectra at finite temperature and density are depicted in figure |l|. 




Figure 1. The photon spectral function R(uj) = —12-k/uj 2 Imfl(a;) at finite temperature 
and p = (left panel) and at T = and finite density (right panel). The qq line in the 
left panel shows the spectral function in the QGP phase for massless quarks. 



4. Dilepton rates 

Eq.(|I|) is now integrated over the space-time evolution of the fireball and the detector 
acceptance of the CERES /NA45 experiment at CERN, using the expressions for ImlT 
from the previous sections. We use simple parametrizations, adapted to the CERES 
conditions, for the volume V(t) and the temperature T(t) J| based on more involved 
calculations utilizing microscopic transport equations. The corresponding rate and the 
data |6| are shown in figure ||. We have considered two cases for the QGP phase: a 
scenario with massive quarks according to eq.([|) and an ideal gas scenario with massless 
quarks. 

The low invariant mass region 0.3 - 0.6 GeV is mainly filled by the strong broadening of 
the p-meson at finite density, contributions from the early QGP phase are subdominant. 
In the quasiparticle approach, the quark masses appear to become heavy close to Tq 
(m q ~ 0.3 — 0.5 GeV), thus there is no radiation from the QGP below the 2m q threshold 
in the invariant mass. 

The uj meson almost loses its quasiparticle structure at high temperatures T > 120 MeV 
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Figure 2. Dilepton rates as a function of the invariant e + e" mass M, calculated in a 
mixed scenario (hadrons below Tq, quarks above Tq = 150 MeV) for 158 AGeV Pb-Au 
collisions at CERES. 



because of the scattering con — > txtx off thermally excited pions. Its peak structure is 
not visible in the dilepton yield and this may act as a signal for medium effects when the 
corresponding experimental resolution is increased. On the contrast, the <fi meson remains 
a distinct peak over the flat QGP spectrum and can be used as both a thermometer and a 
chronometer for the fireball expansion. In the high invariant mass region, the contribution 
from the hadronic phase dies away, and the slope of the dilepton yield is mainly dominated 
by the averaged temperature — (1/T) and is insensitive to the existence of thermal masses. 
We also checked that a possible mixed phase would not leave any distinct traces in the 
dilepton spectrum. 

To uniquely identify signals for the formation of a QGP, it remains a crucial task to 
further investigate temperature and density effects on the hadronic spectra. 
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